Imaging the Space-Time Evolution of High Energy Nucleus-Nucleus Collisions with 

Bremsstrahlung 



J.I. Kapusta* and S.M.H. Wongt 
School of Physics and Astronomy, University of Minnesota, Minneapolis, MN 55455 

(March 2, 1999) 

The bremsstrahlung produced when heavy nuclei collide is estimated for central collisions at the 
Relativistic Heavy Ion Collider. Bremsstrahlung photons with energies below 100 to 200 MeV are 
sufficient to discern the gross features of the space-time evolution of electric charge, if they can be 
separated from other sources of photons experimentally. This is illustrated explicitly by considering 
two very different models, one Bjorken-like, the other Landau-like, both of which are constructed to 
give the same final charge rapidity distribution. 
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I. INTRODUCTION 



The Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory (BNL) will have its first beam 
available for experiments by the end of 1999. A basic quantity of interest is the initial energy density of hot matter 
whose study is the ultimate physics goal. Hadronic measurements may not illuminate the state of the system at 
maximum density for two reasons: First, the system must expand and cool to a sufficiently low density before quarks 
and gluons can bind into independent color singlet hadrons. Second, even after hadronization or confinement phase 
(f-) , transition the hadrons can scatter among themselves, changing their momentum distributions. For these reasons 
hard photons and lepton pairs have been considered better probes of the high temperature plasma of quarks and 
gluons since the production rate of such radiation increases dramatically with temperature. However, high transverse 
momentum photons and lepton pairs can also originate from the lower temperature hadronic phase due to hadron- 
hadron collisions and decays coupled with collective transverse expansion. (Focussing on high mass pairs negates 
this to some extent but eventually the Drell-Yan mechanism will overwhelm these at very high mass.) If it was only 
possible to see these nuclear collisions with the human eye! This is not possible, of course, leaving the next best 
Oh! option: Measurement of the electromagnetic bremsstrahlung produced by the massive deceleration of the nuclear 
pi ■ electric charge. 

CL) , The first theoretical study of nucleus-nucleus bremsstrahlung at relativistic energies, of order 1 GeV per nucleon, 
was carried out by one of the authors fjl]. This was followed by Bjorken and McLerran who were interested in much 
higher beam energies J|], and then by Dumitru, McLerran and Stocker ||. Most recently Jeon, Kapusta, Chikanian 
and Sandweiss studied bremsstrahlung photons with energies less than 5 MeV and a detector design to measure 
them at RHIC. Such an observation would allow for a determination of the final rapidity distribution of the net electric 
charge, which is one measure of transparency and stopping, obviously important for the initial energy density and 
the formation of quark-gluon plasma. It may be argued that a direct observation of all electrically charged hadrons 
would yield the same information, albeit at much higher cost. Nevertheless a cross check would be valuable. 

In this paper, we wish to determine whether it is possible to discern the gross features of the space-time evolution 
of the electric charge using bremsstrahlung photons with energies up to 100 or 200 MeV in the lab frame of RHIC 
(which is also the center-of-mass frame for a collider). Generally a photon of energy 200 MeV would allow one to 
probe distances on the scale of 1 fm and times on the scale of 1 fm/c. At ultrarelativistic beam energies this is not 
quite true due to the extreme forward focussing of the bremsstrahlung. Nevertheless the answer to the question posed 
at the beginning of this paragraph is yes. In section 2, we work out the necessary formulas for bremsstrahlung. In 
section 3, we apply it to two extreme models: the first is Bjorken-like, which has considerable nuclear transparency, 
and the second is Landau-like, which has maximum nuclear stopping. Both models are tuned to produce a flat rapidity 
distribution for the net electric charge so that measurement of hadrons alone could not distinguish them. Section 4 
contains concluding remarks. 
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II. CALCULATION OF BREMSSTRAHLUNG 



Consider a central collision of two equal mass nuclei of charge Z in their center of momentum frame. Denote the 
speed of each beam by vq and the corresponding rapidity by yo, so that the rapidity gap between projectile and target 
nuclei is 2yo- Velocity and rapidity are related by v = tanhy. Before collision each nucleus may be viewed as a 
flattened pancake with negligible thickness since the photons we are considering cannot probe distances smaller than 
1 fm. For example, with a RHIC beam energy of 100 GcV per nucleon the Lorentz 7 = 106 so that the diameter of a 
gold nucleus is contracted to 0.13 fm! 

The formula for computing the classical bremsstrahlung from accelerated charges is well-known ||. For charges 
Qi with coordinates Ti(t), velocities Vi(t) and accelerations &i(t) the intensity and number of photons emitted with 
frequency uj in the direction n are 



d 2 I d 2 N 

LO- 



dujdQ, dcjdQ 



= |A| 2 (1) 



where the amplitude is 
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The sum over charges may of course be replaced by an integral when the charge distribution is viewed as continuous. 



A. Soft bremsstrahlung 

For low frequency photons, the nucleus-nucleus collision appears almost instantaneous. To the extent that the 
transverse rapidities of the outgoing charged particles are small compared to the beam rapidity one may then derive 
the formula j|: 
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Here it is convenient to take n = (sin#, 0, cos 9). The function p represents the charge rapidity distribution at a 
distance r± from the beam axis z and is normalized as: 



dyp(r±,y) = 2, 



(4) 



the 2 arising because the total charge is 2Z. We cannot go further without some knowledge of the distribution p(r±, y). 

It is advantageous to have a simple analytic model with a variable charge rapidity distribution. To this end, suppose 
that p(r±,y) is independent of r±. Then 



where F is a transverse nuclear form factor: 
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where a is the charge per unit area of the nucleus. A solid sphere approximation should be adequate for large nuclei, 
in which case 



Ha) = 



sing 



COS q 



(7) 



where q — ujR sin 9 and R is the nuclear radius. For small angles and low frequencies the nuclear form factor is 
practically equal to one. 

The integral over rapidity can easily be performed for a flat rapidity distribution: p(y) — 1/j/o f° r — j/o < D < Ho- 
Then the photon distribution is: 
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Because vq is close to 1 the distribution deviates strongly from the quadrupole form, peaking near small but nonzero 
9. At RHIC the peak occurs at 9 ~ 1°. The distribution is plotted in Fig. [j] for central gold collisions at RHIC with 
u}« 1/i? = 30 MeV. 
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FIG. 1. The intensity as a function of angle for very soft photons which results from a flat charge rapidity distribution. The 
intensity is independent of frequency for uj <C 1/ FI- 
TS. Hard bremsstrahlung 

All the results above were already reported in |1 and are applicable for very soft photons. Here we will consider 
photons with energies up to 100 to 200 MeV with the expectation that they will reveal the spatial and temporal 
development of the nuclear collision. We will consider only central collisions which can be triggered on experimentally 
with hadronic detectors, such as a calorimeter. 

For this purpose one can certainly invoke specific space-time models of nuclear collisions at RHIC. In this paper we 
will consider simple parameterizations of such collisions which represent extreme limits of such models. The general 
formulas will be worked out here, in this section, and then applied to these models in the next section. 

We will represent a large nucleus, such as gold, by a uniform spherical distribution of nucleons with proper density 
no = 0.155 nucleons/fm 3 . An infinitesimal unit of charge will have a coordinate specified by r± and by the distance I 
parallel to the z-axis as measured from the backside of the nucleus. For given r± the allowed value of I ranges from 
to l ma x{f±) = 2\/R 2 — t\- Summation over all charges i in eq. (2) is then replaced by an integral over and I. 
See Fig. §. 

Each charge element undergoes some deceleration to a final rapidity y. This deceleration will take a finite time 
and can vary from one element to another. In this paper we will fix the final observed charge rapidity distribution 
and vary the space-time evolution to see whether hard bremsstrahlung can provide information not obtainable by 
measuring hadrons only. For definiteness we shall fix the final charge rapidity distribution to be flat. In individual 
nucleon-nucleon collisions it is known to be approximately proportional to coshy in the center-of- momentum frame 
. In order to obtain a flat rapidity distribution in a central nucleus-nucleus collision different charge elements must 
have a rapidity that depends on rj_. The assumed relationship among these variables is 

I = f(rx)y + g(r±_)sinh.y (9) 

where / and g are functions yet to be specified. This means that the last charge element in line will eventually come 
to rest in the center-of-momentum frame. 
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FIG. 2. The charge element coordinate labels in the rest frame of the nucleus. 



The charge per unit area per unit rapidity is: 
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To have a net flat rapidity distribution requires that 
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dr^r ± g(rj_) = 0. 
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As r±_ — » R the distribution should have a minimum at y = 0, in accord with proton-proton scattering, whereas 
when r±_ — > the distribution must have a maximum at y = to produce a distribution that is flat overall. This is 
physically sensible because there ought to be significantly more deceleration in the central core of the nucleus than 
at its outer periphery. A linear extrapolation of nucleon-nucleon scattering to central collisions of gold nuclei at 
RHIC |^||| produces just such a behavior with a net distribution that is approximately flat in rapidity. We choose 
the parameterization g(r±) = gil m ax(r±) — 92^ nax (j' ±) where g\ and <?2 are positive constants. The requirement 
of a net rapidity distribution which is flat results in g2 — 2gi/3i?. The physically sensible requirement that the 
leading charge elements emerge from the collision with the original beam rapidity determines the function f(r±) via 
lmax(r±) = f(r±_)yo + g(r±) sinhyo- Finally, the requirement that the charge rapidity distribution be non-negative 
for all allowable y leads to the condition g\ < 3/(yo cosh yo — sinhyo)- For definiteness we choose the upper limit 
,9i = 3/(yo cosh yo — sinh j/o) since this value brings the rapidity distribution as close to coshy as r±_ — > R as possible, 
consistent with free space nucleon-nucleon scattering. The resulting d 3 Q/d 2 r±dy is plotted in Fig. ||, and is very 
similar to what is obtained in LEXUS (see Fig. 1 of Q). 



For a central collision of equal mass nuclei at RHIC the primary accelerations occur along the direction of the beam 
axis. By azimuthal symmetry the transverse accelerations will, in addition, tend to cancel in the expression for the 
amplitude. The contribution to the amplitude from the projectile will then be: 
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FIG. 3. The charge rapidity distribution, d 3 Q/d 2 r±dy, as it depends on the transverse coordinate r±. 




(12) 



4 



The expression for the contribution from the target nucleus A- will be the same except for the reversal of sign of 
the position, velocity, and acceleration. Since there is a one-to-one relation between the longitudinal label I and the 
rapidity y this amplitude may be written as an integral over y instead of over I as: 

> Ze f R , j / ■ m f V ° j dl(y,r ± ) 



A + =s'm9— — jiq / dr±r± Jq (r±cu sin 9) / dy- 
A Jo Jo 



dy 



dteMi" (t - z(y, t) co ,9)} jr -^— w . (13) 
For a specified model of the deceleration this three-dimensional integral must be evaluated numerically. 



III. TWO EXTREME MODELS 



In this section we will describe two extreme models for the space-time evolution of a central nucleus- nucleus 
collision at RHIC. The real world will undoubtedly be different than either of these. However, our point is that 
hard bremsstrahlung may be used to discern the gross features of the evolving charged matter, something that 
hadron detection alone cannot do because the final charge rapidity distribution in these two models is the same (by 
construction). 



A. A Bjorken-like model 

This model assumes that the nucleons undergo a constant deceleration upon nuclear overlap. It essentially assumes 
transparency of electric charge although the final charge rapidity distribution is assumed to be flat. Because the 
nuclei are so highly Lorentz contracted initially it is an excellent approximation to assume that all charge elements 
are located at z = at the moment t = Q of complete nuclear overlap. A photon with energy less than 200 MeV 
would never be able to distinguish this from a Lorentz contracted pancake of thickness 0.13 fm! For a projectile charge 
element with nuclear coordinates I, r± the position as a function of time is: 

( vot t < 

z(l,r±,t) = I vot + a(l,rj_)t 2 /2 0<t<t f (14) 

[ z(l,r ± ,t f )+v f (t-tf) t f <t. 

The deceleration is 

Vf-Vo tanh y - tanh y 

a{l,r±) = — = • (15) 

*/ */ 

Remember that I and Vf = tanhy are related by eq. (9). The trajectory of a target charge element is obtained 
symmetrically. Note that the charges continue to decelerate even after the nuclei have passed through each other. 
This may be thought of as strings connecting projectile and target charges causing them to decelerate at a distance, 
or by the finite time necessary for particle production. 

The deceleration time tf is allowed to depend on the length of the tube as follows: 

tf(r±) = to + lmax(r±)/co (16) 

where to and cq are constants. Several choices of these parameters will be considered. 



B. A Landau-like model 



This model assumes that the nuclei stop each other within a distance of one Lorentz contracted nuclear diameter, 
then accelerate longitudinally due to the high pressure J8|. The deceleration occurs over such a short time interval 
that photons of 200 MeV or less cannot resolve it. Thereafter our variation of Landau's model assumes that the 
charge elements accelerate uniformly to their final velocities. The trajectories are parameterized as: 

a(l, r±) = ~v S{t) + V -^6{t)6{tf - t) . (17) 

The relationship among I, Vf and y is as above. In a hydrodynamical approach the final time tf is expected to be 
independent of r± but it can be left as a parameter, just as in the Bjorken-like model. 
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C. Comparison of models 



Before comparing the bremsstrahhmg produced in these two models let us compare the space-time evolutions. In 
Fig. ^ we plot cells bounded by lines of constant I and r± in the x — z plane as a function of time. Each cell contains 
the same amount of charge. The left panel shows the Bjorken-like model and the right panel shows the Landau- 
like model. Both models have tf = 1 fm/c; that is, the acceleration time is independent of r±. There are several 
noteworthy features: The Bjorken-like model always has a gap in the beam direction because it takes a finite time for 
the charges to decelerate. The charges that finally come to rest in this frame of reference require 1 fm/c to decelerate 
from v = ±vq to v — and by then they have traveled 0.5 fm. In contrast, the charges in the Landau-like model stop 
instantaneously and then accelerate to their final velocities. Therefore there is no gap in the beam direction. Note that 
the shaded regions indicate where the net charge is nonzero. The gap in the central region in the Bjorken-like model 
does not signify empty space but rather a region where there is zero net charge; it may be occupied by quark-gluon 
plasma or hadronic matter or strings or something else. 



Bjorken-like Landau-like 
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FIG. 4. Space-time evolution in a Bjorken-like model compared to that of a Landau-like model both with tf = 1 fm/c. The 
curves indicate contours of constant I and r± and so each cell, as delineated this way, contains the same amount of electric 
charge. 

In Fig. U we show the space-time evolution for to = 0.5 fm/c and Co = 1 in the Bjorken-like model and for 
to = 10 fm/c and cq — oo in the Landau-like model. We chose the acceleration time in the Landau-like model to be 
independent of rj_ because this seems more relevant for a hydrodynamical model. In the Bjorken-like model we chose 
the minimum time to = 0.5 fm/c because, even in free space, nucleon-nucleon collisions must occur over some finite 
time interval. The value of 0.5 fm/c is suggested by the analysis of Gale, Jeon and Kapusta |9| of Drell-Yan and J/ijj 
production in very high energy proton-nucleus collisions. The central cores of the gold nuclei then require 13.9 fm/c 
to decelerate to their final velocities. Hence the average deceleration time ought to be comparable to the 10 fm/c 
chosen for the Landau-like model. Apart from the quantitative difference in the space-time evolution, the gap in the 



G 



Bjorken-like model is now more like an ellipse than a slab. 



Bjorken-like 




Landau-like 




FIG. 5. Same as Fig. |4| except that the Bjorken-like model has at/ varying from 0.5 fm/c at rj 
and the Landau-like model has fixed tf = 10 fm/c. 



R to 13.9 fm/c at r x = 0, 



The intensity distribution in the Bjorken-like model is plotted in Fig. ^ for the angles 1°, 5°, 10° and 20°. Three 
values of cq have been chosen for display: 0.5, 1 and oo. The value of to has been fixed at 0.5 fm/c in all cases, 
although the results are practically the same with any to < 1 fm/c. The three values of Co correspond to final times 
for the central cores of 27.3, 13.9 and 0.5 fm/c, respectively. (Be careful to note the false zeroes in three of the four 
panels.) The dependence on Co is very weak considering the wide range of resulting tf. This may be understood as a 
relativistic effect. Because of the relation v = tanh y, most of the outgoing charges still move with a speed nearly equal 
to that of light. Therefore the phase in eqs. (12-13) is approximately iu>t(l — cos 8), and for it to become comparable 
to 1 requires an lu much bigger than naively expected. The separation of the curves increases with 8 following this 
phase factor. The primary falloff of the intensity with increasing photon energy is caused by the transverse form 
factor described by eq. (7). The falloff increases with angle because the argument of the form factor is q = ojRsmd. 

The intensity distribution in the Landau-like model is plotted in Fig. fj] for the same set of angles. Four illustrative 
values of tf have been chosen: 1, 2, 10 and 20 fm/c. In this model there is a striking dependence on tf. The oscil- 
latory behavior for large tf is due to the interference between the initial instantaneous stopping with the subsequent 
acceleration to the final velocities. At 10° and with tf = 20 fm/c, for example, the enhancement due to constructive 
interference is an order of magnitude at lu — 60 MeV as compared to the soft photon limit. 
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Landau-like 
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FIG. 7. The intensity distribution in the Landau-like model as a function of photon energy for various angles. The four 
curves correspond to different acceleration times. 
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The oscillatory behavior which arises in the Landau-like model but not in the Bjorken-like model is easy to un- 
derstand in a simplified version of these models. Let us shrink each nucleus to a point charge. In the Bjorken-like 
version, assume that the charges decelerate from ±vq to ±Vf instantaneously. The intensity distribution is: 

ckj<m~ n 2 sm (i -v$ cos 2 0) 2 (i -v* cos 2 ey ! (18) 



which has no frequency dependence. In the Landau-like version assume that the charges stop instantaneously remain 
at rest for a time r, then instantly accelerate to final velocities ±v/. This intensity distribution is: 

d 2 I Z 2 a 9 9 \ [vq - V 2 f cos ujt - v 2 v 2 f (l - cos wr) cos 2 Of 
■ sin 6 cos ' 



dudVL tt 2 1 (1 - w 2 cos 2 (9) 2 (l -v 2 cos 2 6) 2 

(19) 
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vj sm ujt 
(1 — v 2 cos 2 9 



The period of oscillation is just the time delay between stopping and sudden acceleration. This is essentially what 
happens in the Landau-like model, modified by the smearing which follows from a continuous acceleration rather than 
a time delayed one, and by the relativistic effect in the phase factor as mentioned before. 



IV. CONCLUSION 



In this paper we have studied two very different and extreme models of the space-time evolution of the charges in 
central nucleus-nucleus collisions at a RHIC energy of 100 GeV per nucleon in the center of mass frame. The models 
were constructed to yield the same final rapidity distribution of the net outgoing electric charge. This means that 
measurement of hadron spectra alone cannot discern them. Very low energy bremsstrahlung cannot discern them 
either because it depends only on the difference between the incoming and outgoing currents, which are identical. 
Bremsstrahlung photons of higher energy can distinguish these extreme models. For the Landau-like model photons 
of energy 100 MeV or so is quite sufficient to infer the time scale for expansion. Bremsstrahlung photons with energy 
up to 200 MeV may be able to infer the deceleration time in Bjorken-like models but one needs to look at angles 
much greater than where the peak of the spectrum occurs, which is 1 degree. Whether it is really possible to design 
a detector that can separate other sources of photons of this energy is an open question, but at least we have shown 
that the physics to be learned is interesting, perhaps justifying an experimental effort [fl0|| . 
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